Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA     

Quansheng Liu  Weiqin Wang  Hao Ma 《国际地质力学数值与分析法杂志》2020,44(2):208-238

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.  相似文献   

重载铁路水泥改良膨胀土路基填料动弹模量及阻尼比研究     

商拥辉  徐林荣  黄亚黎  陈钊锋 《工程地质学报》2020,28(1):103-110

动弹模量与阻尼比是土动力学分析中的重要力学参数,考虑重载铁路荷载特征定量分析水泥改良膨胀土的动模量和阻尼比的较少。依托蒙西至华中地区铁路煤运通道(简称蒙-华铁路)工程为背景,采用南阳邓州市大山寨膨胀土,通过在不同频率、围压、固结比及动应力幅值下的持续振动三轴试验,研究了水泥掺量3%和5%水泥改良膨胀土的动弹模量及阻尼比,并与膨胀土素土进行对比分析。结果表明:水泥掺量3%和5%改良膨胀土的最大动弹模量约为膨胀土素土的3~4倍;在动弹模量-应变曲线中,动应变小于0. 002时表现为陡降段,动弹模量随动应变增长降幅达70%,而动应变大于0. 002时降幅较小,动弹模量随动应变增长趋于稳定;动弹模量随围压、频率、水泥掺量增加而增大,阻尼比随围压、固结比增加而减小;低应变水平下,固结比与动模量成正相关关系,高应变水平下,固结比与动弹模量成负相关关系。同时,对动弹模量及阻尼比进行了归一化分析,建立了估算动弹模量及阻尼比的经验公式。  相似文献   

Study of irregular wave run-up over fringing reefs based on a shock-capturing Boussinesq model     

《Applied Ocean Research》2019

This paper presents the results of a parametric study of irregular wave run-up over fringing reefs using the shock-capturing Boussinesq wave model Funwave-TVD to better understand the role of fringing reefs in the mitigation of wave-driven flooding. Laboratory experiments were newly performed with a typical fringing reef profile and typical hydrodynamic conditions to validate the model. Experimental data shows irregular wave run-ups are dominated by the low-frequency motions and confirms the run-up resonant phenomenon over the back-reef slope, which has been revealed in previous numerical studies. It is demonstrated that irregular wave evolution and run-up over fringing reefs are reasonably reproduced by the present model with a proper grid size. However, the infragravity run-up height and highest 2% run-up height over the back-reef slope are under-predicted due to the underestimation of the infragravity wave height over the reef flat. The validated model was then utilized to model irregular wave transformations and run-ups under different conditions. Through a series of numerical experiments, the effects of key hydrodynamic and reef geometry parameters, including the reef flat width, water depth over the reef flat, fore-reef slope angle and back-reef slope angle, on the irregular wave run-up were investigated. Variations of spectral components of irregular wave run-ups were examined to better understand the physical process underlying the effect of each parameter.  相似文献   

GNSS多天线基线网单历元模糊度同步解算法     

汤佳明  柴艳菊  闻德保  欧吉坤 《大地测量与地球动力学》2019,39(3):262-268

提出多天线基线网单历元模糊度同步解算法，其具体实现步骤为：1）根据宽巷组合模糊度易于固定的优点，采用附加已知基线长度约束法同步解算各基线的宽巷模糊度，得到dm~cm级精度的近似基线分量；2）将解算得到的各近似基线分量作为约束，同步解算各基线的基频模糊度，以获取mm级精度的基线矢量。该方法的关键在于检验各历元宽巷模糊度解算的正确性，以获取可靠的近似基线分量，为解算各基线的基频模糊度提供准确的基线先验信息。由于动态情况下各历元观测信息比较少，单纯依赖ratio检验不可靠，提出结合基线误差、单位权中误差、基线网模糊度闭合环及ratio值等对多组宽巷模糊度进行检验，避免ratio值设置不当导致模糊度检验中发生纳伪和弃真问题。实测数据结果证明，该处理方法使得模糊度解算的成功率提高1%~2%，可以获取移动平台更丰富的导航信息，提高其服务能力。  相似文献   

唐山地区波速比变化特征研究     

王亚茹  宋昭  刘晓丹  郭蕾  王妍  王时  宫猛  王想 《大地测量与地球动力学》2019,39(7):686-690

选取唐山地区2008~2018年震相数据，利用单台多震和达法和多台多震和达法分别计算波速比，结合研究区内的地震活动对波速比的变化特征进行研究。结果显示，多台多震和达法得到的波速比结果较为稳定，而单台多震和达法得到的结果变化幅度大，显示更多细节；唐山地区ML≥4.5地震发生前单台波速比存在不同程度的异常，异常台站的方位与地震具有一定的对应性。  相似文献   

月球辐射状撞击坑经度方向分布特征与同步旋转轨道状态模拟分析          下载免费PDF全文

邵先远  李斐  鄢建国  叶茂  郝卫峰 《地球物理学报》2017,60(4):1321-1331

月球目前的同步旋转轨道状态使得形成的撞击坑分布满足一定的不对称性.本文利用最新的LRO影像和地形数据,结合早期的Clementine影响数据,分析了月球辐射状撞击坑经度方向成坑率分布,结果表明所识别的辐射状撞击坑的.年龄为O.9 Ga内,西一东半球方向存在明显不对称性,比值约为1.35~1.53;同时利用辐射状撞击坑模拟了0.9 Ga内的同步旋转轨道的平均状态,得到该时间内月球的向点-背点为70°W—110°E附近;最后利用撞击坑数据库资料对月球形成以来各地质世纪时间尺度内的撞击坑进行了成坑率分布分析,结论表明月球在大爆炸中后期间可能处于过近似的同步轨道旋转状态,但是与现在的轨道运行状态相反,月球在之后的地质时期内经历了翻转,其诱因可能是雨海和东方海盆地遭受的撞击.  相似文献   

基于动态图像变形的PP与PS波层位直接匹配   总被引：1，自引：1，他引：0       下载免费PDF全文

蒋雪珍  芦俊  王赟 《地球物理学报》2017,60(3):1106-1117

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分:首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.  相似文献   

CAP1400核岛结构缩尺模型对地震反应影响分析     

李小军  韩杰  王晓辉  贺秋梅 《震灾防御技术》2017,12(2):241-253

针对核电厂CAP1400核岛结构地震反应问题，构建了核岛屏蔽厂房和辅助厂房整体结构的3个分析模型:原型和1/16、1/40缩尺模型，并在AP000谱和RG1.60谱地震动输入下进行了有限元模拟对比分析，探讨了振动台试验模型缩尺处理的合理性和精确性。研究表明，基于缩尺模型得到的结构自振频率相对于原型结构模型有所降低，降低幅度在8.5%以内；结构模型的缩尺对结构反应峰值加速度和高频（大于3Hz）加速度反应谱的影响较为显著，但对较低频（小于3Hz）的加速度反应谱影响较小；模型缩尺对结构不同方向反应的影响中，刚度越大的方向其影响越大。进一步将结构模型数值模拟结果与1/16缩尺模型的振动台试验结果进行了比较分析，试验给出的结构自振频率远低于模型数值模拟结果，但原型和1/16缩尺模型数值模拟得到的结构反应均与试验结果较为接近。基于模型数值模拟和振动台试验研究，认为对于缩尺比1/16或更大的模型可以忽略模型的缩尺效应。  相似文献   

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion          下载免费PDF全文

Jiazhong Qian  Lei Ma  Hongbin Zhan  Qiankun Luo  Xiao Wang  Mu Wang 《水文研究》2016,30(11):1718-1726

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Helmert方差分量估计在GPS/GLONASS/BDS组合定位权比确定中的应用     

刘金海  涂锐  张睿  张鹏飞  卢晓春 《大地测量与地球动力学》2018,38(6):568-570

因GNSS系统间观测噪声、轨道精度的差异，采用经验权比进行组合定位难以得到最优结果。基于此，在GPS/GLONASS/BDS组合定位中引入Helmert方差分量估计，对GPS/GLONASS/BDS组合单点定位和基线解算中各系统观测值进行合理定权。实验表明，采用该方法确定的伪距观测值最佳权比为5∶1∶1，相位观测值最佳权比为1∶1∶1，有效提高了GPS/GLONASS/BDS组合定位的精度和可靠性。  相似文献